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ABSTRACT 
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We discuss Green's-function solutions of the equation for a geometrically thin, 

axisymmetric Keplerian accretion dis c with a viscosity prescr i ption y oc R n . The 
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mathematical problem was solved by iLynden-Bell fc Pringld (119741 ) for the spe- 
r^ cial cases with boundary conditions of zero viscous torque and zero mass flow 

Y*\ at the disc center. While it has been widely established that the observational 

appearance of astrophysical discs depend on the physical size of the central ob- 
ject (s), exact time-dependent solutions with boundary conditions imposed at 
finite radius have not been published for a general value of the power-law index 
n. We derive exact Green's-function solutions that satisfy either a zero-torque or 
a zero-flux condition at a nonzero inner boundary R\ n > 0, for an arbitrary initial 
j — | surface density profile. Whereas the viscously dissipated power diverges at the 

disc center for the previously known solutions with i? in = 0, the new solutions 
with R in > have finite expressions for the disc luminosity that agree, in the 



O 

limit t — y oo, with standard expressions for steady-state disc luminosities. The 

new solutions are applicable to the evolution of the innermost regions of thin 
accretion discs. 



Subject headings: accretion, accretion discs 



1. Introduction 



Since its emergence in t he 1970's ( iShakura &: Sunyaevlll973t iNovikov &: Thorndll973 



Lynden-Bell fc Pringld 119741 ) . the theory of astrophysical accretion discs has been applied 
to explain the emission properties of active galactic nuclei, X-ray binaries, cataclysmic bi- 
naries, supernovae, gamma-ray bursts, and the electromagnetic signatures of mergers of 
supermassive black holes; to study planetary and star formation; and to model the evolution 
of binary and planetary systems. Because accretion discs onto compact objects can dissipate 
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much larger fractions of baryonic rest-mass energies than nuclear reactions, they are often 
associated with some of the most energetic astrophysical processes in the universe. 

If the local gravitational potential is dominated by a central compact object or a compact 
binary, and if the timescale for the viscous dissipation of energy is longer than the orbital 
timescale, then the accretion flow near the center of the potential is expected to be nearly 
axisymmetric. If the gas is able to cool efficiently, then the flow will also be geometrically 
thin, and one only needs the radial coordinate to describe the mass distribution in the disc 
(any relevant vertical struct ure can be integrated or ave raged over the disc height). The 
partial differential equation f lLynden-Bell fc Pringldll974l . henceforth LP74) 



d_ 
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E{R,t) 



RdR 



Rl/2 M ^ R ^ 



is obtained by combining the equations of mass conservation and angular momentum, and 
describes the surface density evolution of a thin Keplerian accretion disc due to kinematic 
viscosity v. 

In general, the viscosity v depends on the surface density E and equation (pQ) is nonlinear. 
If, however, v is only a function of radius, then the equation is linear and much more amenable 
to analytic methods. In particular, a solution that makes use of a Green's function G, 



/■oo 

E(R, t)= G(R, R', t) E(R', t = 0) dfl, 

jR in 



(2) 



gives the solution E for any t > given an arbitrary profile H(R,t = 0) and an inner 
boundary condition imposed at R{ n . A distinct advantage of the formalism is that it gives 
the solution T,(R,t) through a single ordinary integral, whereas a finite-difference algorithm 
would require the computation of the profile at intermediate times. Another advantage is 
that the initial density profile need not be differentiable. 



For a power-law viscosity v oc R n , iLiistl ( 119521 ) and ILP74I derived analytic Green's 
functions that satisfy a boundary condition of either zero-torque or zero-mass-flux at the 
coordinate origin, i.e., for the case R in = 0. In reality, however, the objects at the center 
of astrophysical accretion discs have a finite size to which the observational appearance of 
the disc is sensitive: e.g., the luminosity, spectral hardness, and variability timescales of 
black hole discs depend strongly on the radius of the innermost stable orbit, and those of 
circumbinary discs depend on where the inner disc is truncated by the central tidal torques. 
Green's functions with R in = do not capture the time-dependent behavior of accretion 
discs close to the central object. In fact, in solutions obtained with such Green's functions, 
the integral for the total power viscously dissipated in the center of the disc diverges. 
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Despite the astrophysical relevance of Green's functions to the thin accretion disc equa- 
tion with boundary con ditions imposed at a finite radius, such solutions have not been 
published. iPringld (Il99lf ) derived the Green's function with a zero-flux boundary condition 
at a nonzero radius in the special case n — 1, and noted the "extreme algebraic complexity" 
involved in calculating a more general solution with R m > 0. Time-dependent models of 
accretion flows have continued t o employ solutions t hat correspond to the c entral objects 
having zero physical size (e.g., iMetzger et al.l l2008t iTanaka fc Menoul 120101 1 . In order to 
calculate a convergent disc luminosity and spectrum, such models typically approximate an- 
alytically the effects of an inner boundary condition, e.g., by truncating the disc profile at 
an artificially imposed radius. 

In this paper we derive exact Green's functions for equation ([T]) for boundary conditions 
imposed at a finite radius, for any power-law viscosity v oc R n with n < 2. We show 
that mathematical difficulties can be minimized with the aid of the ap propriate integral 
transform techniques, namely the Weber transform (|Titchmarshl 119231 ) and the recently 
proved generalized Weber transform (IZhang fc Tongll2007l ). We present two specific solutions 
of astrophysical interest: the solution with zero torque at a radius R m > 0, which is of interest 
for accretion discs around black holes and slowly rotating stars; and the solution with zero 
mass flow at R in > 0, which is applicable to accretion flows that accumulate mass at the 
disc center due to the injection of angular momentum from the tidal torques of a binary or 
perhaps the strong magnetic field of the central object. 

Thi s paper is o rgani zed as follows. In §2, we review the Green's function solutions, de- 
rived by lLiistl ( 119521 ) and lLP74l . for the thin-disc equation with boundary conditions imposed 
at the origin. In §3, we derive the new Green's function solutions, which impose boundary 
conditions at a finite inner boundary radius. We offer our conclusions in §4. 



2. Green's-Function Solutions with Boundary Conditions at R = 



In the special case where the viscosity is a radial power law, v oc R n , and assuming a 
separable ansatz of the form E(i?, t) = R p a(R) exp(— At), where p and A are real numbers 
and a is an arbitrary function of R, equation ([[]) can be rewritten as the Bessel differential 
equation: 



R : 



d 2 a 
dR 2 



2p + 2n + 
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0. 



(3) 
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Above, s = vR n is a constant. With the choices p = n — 1/4 and A = 3sk 2 , equation ([3]) 
has the general solution 

a k (R) = R- 2n [A(k)J t (ky) + B(k)Y E (ky)} . (4) 

Above, k is an arbitrary mode of the solution; A(k) and B(k) are the mode weights; £ = 
(4 - 2ra)" 1 > 0; y(R) = R^~ n ^/(1 - n/2); and J^ and Y e are the Bessel functions of the 
first and second kinds, respectively, and of order t. If I is not an integer, then Yg above may 
be replaced without loss of generality by J_£. Integrating the fundamental solution across 
all possible fc-modes gives the solution: 

poo 

E{R,t)= R- n - 1/4 [A(k)J e (ky) + B(k)Y E (ky)} exp(-3sk 2 t) dk. (5) 

Jo 

The mode-weighting functions A(k) and B(k) are determined by the boundary condi- 
tions and the initial surface density profile H(R,t = 0). Our goal is to rewrite equation (jSJ) 
in the Green's function form (equation [2]) and to write down an explicit symbolic expression 
for the Green's function G(R,R',i). Throughout this paper, we will employ the following 

strategy: 

1. Using the boundary condition, find an analytic relationship between the mode weights 
A(k) and B(k). 

2. Identify the appropriate integral transform to express the mode weights in terms of 
the initial profile T,(R,t = 0). 

3. Insert the time-dependence exp(— 3sk 2 t) and integrate over all modes to find the 
Green's function. 

4. Derive analytic expressions for the asymptotic disc behavior at late times and small 
radii. 

Before deriving the solut ions w ith boundary conditions at finite radius, we begin by reviewing 



the Green's functions of ILP74J with boundary conditions at the coordinate origin. 



2.1. Zero torque at R m = 

An inner boundary condition with zero central torque is of astrophysical interest as it 
can be used to describe accretion onto a black hole or a slowly rotating star, at radii much 
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larger than the radius of innermost circular orbit or the stellar surface, respectively. The 
radial torque density g in the disc due to viscous shear is 

g(R,t)^^R^oc^R^, (6) 

where Qk is the Keplerian angular velocity of the orbit. 

Because the functions Ji and Yi have the asymptotic behaviors Ji(ky) oc y e oc i? 1//4 and 
Ye(ky) oc y~ e oc i?~ 1//4 near the origin, at small radii the mode weight A(k) will contribute 
to the behavior g oc R 1 ^ 2 while B(k) will contribute to g — constant. Therefore, for the 
solution to have zero viscous torque at R = the function B(k) must be identically zero. 

We may relate the surface density distribution at t — and the weight A(k) via the 
integral equation 

/"OO 

E(R,t = 0) = R-"- 1 /* / A(k)J £ (ky) dk, (7) 

Jo 



which may be solved with the use of the Hankel integral transform (e.g., IOgilvidl2005l ). 
A Hankel transform pair of order £ satisfies 

poo 

<j) e (x) = / $*(fc) Ji{kx) k dk, (8) 

Jo 

oo 



$i(k) = / 4>e(x) Je(kx) x dx. (9) 

Jo 

For the problem at hand, the suitable transform pair is 

R n+1 / 4 j:(R,t = 0) = / [Aityk- 1 ] J e {ky) k dk, (10) 

Jo 

/"OO 

A(k)k~ l = / [i2 n+1/4 £(i2,£ = 0)] J e (ky) y dy. (11) 

Jo 

Combining them gives us A(k): 

A(k)=(l-iy 1 J e °2tf,0) Hky')kR*/UR' (12) 

Inserting equation (Tl2|) and B{k) = into equation (j5]), we obtain 

E(i2, t) = 1 1 - - ) i?- n -V4 / ^'5/4 / E ( R ,^ t = ) j^^ j^ ky j exp (_ 3sfc 2^ fc dfc dR/ 

(13) 
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To pose the solution in terms of a Green's function G(R,R',t) (equation [2]), we write 



G(R, R , t) 



n 



1 r°° 

--) ^-"-1/^5/4 ^ j^/) j^^ exp(-3sA; 2 t) jfe rfA; 



,R~ 9 / 4 R' 5 / 4 

^- n) ^iRT h 



2 (#/#) 



l-re/2' 



T(R) 



exp 



1 + CgAg) 

t(R) 



2-ra 



(14) 



Above, if is the modified Bessel function of the first kind, and we have substituted t(R) = 
12(1 -n/2) 2 R n ~ 2 st = 8(l-n/2) 2 [t/t u (R)}, where t„{R) = (2/3)R 2 /u(R) is the local viscous 
timescale at R. 

Although the Green's function allows for the calculation of £(/?,£) for arbitrary initial 
surface density profiles, it is instructive to study the case where the initial surface density is 
a Dirac 5 function, 

E{R,t = 0) = S S(R-Ro) Ro, (15) 

for which the solution is (by definition) the Green's function itself. The integral over ra- 
dius in equation (j2]) becomes trivial and many behaviors of the solution may be expressed 
analytically. Because any initial surface density profile can be described as a superposition 
of 5-functions, studying this special case will help illuminate the general behavior of all 
solutions. 

We may evaluate the asymptotic behavior at late times and small radii by noting that 
for small argument z £ 0.2 VT+1, h(z) « {z/2) l /T(l + 1). We find 



iy,0 



where t ufl = t u (R Q ). 

Thus, for these solutions the inward radial mass flow, 

d 
M{R) = -2'kRT.vr = 6nR 1/2 — (z/Si? 1/2 ) 

becomes radially constant near the origin and at late times: 



(16) 



;i7) 



2 — n 
M{t>U{R))^—^--M 

Above, we have defined M = 37rz/(i? )E . 



2 / t v n 
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The power per unit area that is locally viscously dissipated from each face of the disc 
is F = (9/8)^Ef2 2 . The total power dissipated near the center of the disc diverges: 



f R 9tt f R 

L(R,t> t„(R)) = / ^-v(R')Z(R' , t)n 2 (R')R' dR' oc / R'~ 2 dR' 
Jo 2 J 



(19) 
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Astrophysical accretion flows do not extend to zero radius, and thus in practice one may 
truncate the disc at some plausible boundary radius, for example the radius of innermost 
stable circular or bit for a disc around a black hole, by approximating the effects of a finite 
boundary radius ( 1LP74I ). 



Although we have used a 5 function for demonstrative purposes, the quantities Ro 
and E that set the physical scale and normalization of the initial surface density profile, 
respectively, are arbitrary. The asymptotic behaviors noted above hold for any initial surface 
density profile: at late times, the surface density profile approaches £ oc R~ n , M becomes 
radially constant, and the disc luminosity L formally diverges at the center. 



At early times and large radii, such that t <C ^t u (R) t u (R), we may use the fact that 
U(z 3> 1) ~ exp(x)/v / 2vri to find 



G (t < y/t u (R) t^Rffj 



1 



V^iR) 



cxp 



1 - (R'/R) 



l-ra/2 



-i 2 



T(R) 



R >\ (3/4)(l+n) 
Rj 



d 

~dR' 



r_ 

~R 
(20) 



In Figure [U we plot the solution £(/?,£) and the radial mass flow M(R,t), for the 5- 
function initial condition (equation IT5l) . and for viscosity power-law index values n — 0.1 
and n — 1. In both cases, we see the power- law behavior from equation f lT6|) near the origin 
as the solution approaches t ~ t v <Q. The disc spreads as the gas at inner annuli loses angular 
momentum to the gas at outer annuli. The gas initially accumulates near the origin, then 
becomes diffuse as mass is lost into the origin. 



fl-n/2" 



2.2. Zero mass flow at R„ 







If the accretion flow has a sufficiently strong central source of angular momentum, then 
the gas will be unable to flow in, and instead accumulate near th e origin. Such solutions can 
be used to describe astrophysical discs around a compact binary (IPringlell99ll ). and perhaps 
those arou nd compact ob jects with strong central magnetic fields (1LP74I ). For circumbinary 
thin discs, iPringld (Il99l[ ) demonstrated that such a boundary condition characterizes quite 
well the effects of an explicit central torque term. 

In general, the mass flow has the behavior 

M oc / R 1/2 —[A(k)J e (ky)R 1/A + B(k)Yi(ky)R 1/A ]exp(-3sk 2 t)dk. (21) 

Jo "R 

We have seen above that for solutions with B(k) = the mass flow is radially constant 



and finite near the origin. On the other hand, because Yi(ky) oc i? -1 / 4 near the origin, the 
weights B{k) will all contribute no mass flow there; so for zero mass flow at R in = 0, we 
require A(k) = 0. 

We note that because the surface density will have a power-law £ oc Rr l l 2 ~ n at the 
origin, for the mass contained in the disc to converge n must be less than 3/2. Thus, for 
physically realistic solutions with zero mass flow at the origin, i cannot be an integer. It 
follows that in this case Yt in equation fl5]) may be replaced by J_^ without loss of generality. 
Then the Green's function for this case is derived in exactly the same fashion as in the 
previous case, the only difference being that the order of the Hankel transforms has the 
opposite sign. We obtain: 



D-9/4D/5/4 

G(R,R',t) = (2-n) n " U 

t{R) 



2 (R'/R) 



l-n/2 



r(R) 



exp 



1 + (R'/R) 
t(R) 



2-n 



(22) 



As before, we evaluate the late-time behavior for the (^-function initial condition (equa- 
tion [T5l) at small radii: 



E(R,t>t u (R)) 



n 



Enl^ 



-71.-1/2 



1-™ 

2 



2 t 
tv,0 



-i+e 



(23) 



T(l-t) u V#o, 

From the above expression it is clear that the boundary condition is satisfied: M oc d(u'ER 1 ^ 2 ) — > 
in the limit R — > 0. Just as we found for the zero-torque boundary condition, the formal 
expression for the power dissipated at the disc center diverges for the zero-flux solution, with 
L(R <Ro,t> U, ) cc j Ro R-^l 2 dR. 

The asymptotic behavior at early times and large radii is unaffected by the order of 
the function Ii\ it is described by equation [201 Indeed, the inner boundary condition should 
have no effect on the disc at large radii. 

Figure [2] shows the evolution of the surface density and the radial mass flow for the 
boundary condition M(R = 0) = 0. At early times, the behavior is nearly identical to the 
zero-torque boundary case. At late times, the zero-flux boundary condition causes the gas 
to accumulate instead of being lost to the origin. The central mass concentration reaches a 
maximum, then decreases as the disc begins to spread outward. 



3. Green's-Function Solutions with Boundary Conditions at Finite Radii 



As we have seen above, Green's-function solutions of thin accretion discs with i?; n = 
have divergent expressions for the dissipated power, and thus the innermost surface den- 
sity profile must be manipulated to obtain physically realistic disc luminosities. Analytic 
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treatm ent of the case with finite boundary radius was briefly discussed in ILP74I and iPringk 
( 1X99 lh . but to the author's knowledge explicit solutions have never before been published. 
We show below that the Green's functions for finite boundary radii can be derived with the 
aid of the appropriate integral transform techniques, and that they can be represented as 
ordinary integrals of analytic functions. 



3.1. Zero Torque at R m > 

We wish to solve the problem as in §2.11 but with R m > 0, i.e. 



g(R iQ ) oc EuR 1/2 oc Z(R in )R 



71+1/2 



0. 



R — -fxin 



(24) 



We may relate the mode weights A(k) and B{k) by requiring that every mode of the 
solution satisfy the boundary condition, i.e.: 

A(k)J e (ky m ) + B(k)Y e (ky in ) = 0, (25) 

where y m = y(R in ). Substituting C{k) = A(k)/Yi(ky in ) = —B(k)/Ji(ky in ), we obtain 

/•oo 

E(R,t) = / C{k)R- n - 1/4 [J i {ky)Y i {ky in )-Y i {ky)J i {ky in )}exp{-3sk 2 t) dk. (26) 

Jo 



The function C(k) may be evaluated with the use of the Weber integral transform 
(JTitchmarshlll923l ). A Weber transform pair satisfies 



f°° J e (Kx)Y e (K) - Y e (Kx)J e (K) 

Mx) = Jo M) Jj{n) + Yl{n) KdK ' 



$/(k) 



4>e(x) [Ji(kx)Yi(k) — Ye(Kx)Je(K)] x dx. 



(27) 
(28) 



Proceeding as before, we construct the pair 

i? n+1/4 S( J R,t = 0)= / 
Jo 

POD 

C^k- 1 = / [R n+l/A E(R,t = 0)] [Ji(kx)Yi(k) - Yi(kx)Ji(k)] x dx. 



f°° r r( v -11 Je(Kx)Y e (K)-Y t (Kx)J t (K) 

[ {) J jM+WW) 



(29) 
(30) 



Above, we have substituted x = y/y- m > 1 and k = ky in . Note the lower limit of 
integration in equation (130]) is nonzero to account for the finite boundary radius. Combin- 
ing equations ( )29l and ( 1301) to eliminate C(k), and inserting the time- dependence factor 
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exp(— 3sk t) = exp[— 2(1 — n/2) k t/t U: - m ] where t u ^ n = t u (R- m ), we obtain our new Green's 
function: 

G{R, R , £) = 



n 



R-n-l/4 Rf 5/A R n-2 

[J e (Kx)Y e (K) - Y t (Kx)Jt(K)] [Je(Kx')Y e (K) - Y^kx 1 ) J £ (k)} 



x exp 



1- - 
2 



1\2 2 t 



t„ 



K d,K. 



(31) 



Whereas the integral over k in equation ( TH| has an analytic solution, to the author's knowl- 
edge there is no analytic expression for the integral in equation ( |3T1) . Nonetheless, equation 
( 13TJ) gives an exact expression for the Green's function. While it is somewhat more unwieldy 
than the solutions for R in = 0, the additional computational cost of an ordinary integral is 
not likely to be a significant practical barrier, e.g. one could tabulate the integral in terms of 
the quantities x, x' and tjt v ; xn . The boundary condition has little effect at large radii, so in 
practice the behavior far from the boundary is well approximated by the R m = solutions. 

The Green's function in equation (I14p do es have a close d-form expression for the special 
case n — 1 (i.e., £ = 1/2). As noted by iPringld (|199ll ). in this case the Bessel func- 
tions become easier to handle analytically, with J\/2{x) = a/7t/2 a^^sinx and Yi/2(x) = 
— \/ixj2 x^ 1 / 2 cosx. For this value of n we obtain for our Green's function 



G(R,R',t) 



7T 



Rr 



R 



5/4 



[X X 



R-^ 2 R'R; 



-1/2 



2V27T 



-1/2 



exp 



sin [k(x — 1)] sin [n{x — 1)] exp 



K 



(x — x 1 ) t 



exp 



[X + X 



2 t ■ 
2) 2 t 



(Ik 



(32) 



Note that in the case n = 1, x and x' are simply y/R/R- m and ^JR' /R m , respectively. 

For general values of n, we can evaluate the behavior at late times t > £„ )0 > t u ^ n by 
noting that in this regime only the modes k 2 ^ 1 contribute to the integral in equation (13~T|) . 
For the central region R ^ Rq at late times, we obtain the following analytic expression for 
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the ^-function initial condition: 

1 — n 
E(R,t>U(R)) 




We see that the Green's function explicitly gives the asymptotic behavior £ oc R "(1 
y/R- m /R), which has be en used extensively for s olutions of accretion discs near zero-torque 



boundary surfaces (e.g.. lLP74l . lFrank et al.ll2002l )rl This behavior near the boundary and at 



late times is general for any initial surface density profile; it is insensitive to the values of 
S and Ro, and arises for any nonzero R[ n . This qualitative difference in the inner disc from 
the i? in = case also gives a convergent value for the power dissipated in the central disc. 
We obtain: 



L(R<R t ,t>t u (R)) 



GMM S! 

2-Rin 



, Rji 

'lit 



Rjn 
Rt 



3/2 



(34) 



where M ss is the mass flow quantity in equation ( ITS1) . and Rt is the radius where t = t„(R), 
inside which the disc has had sufficient time to approach the asymptotic solution. In the 
limit Ro 3> R m and t ^> t u0 , M ss may be interpreted as the mass supply rate into the center 
of the disc from arbitrarily large radii. In this limit the above expression agrees precisely 
with the standard expression for the luminosity of a steady-state thin accretion disc. 

We show in Figure [3] the exact solutions for the <5-function initial condition, with the no- 
tor que bo undary condition imposed at _R in = Rq/5. The qualitative evolution is as predicted 
by ILP74J : at early times, far from the boundary, the disc spreads inward in very much the 
same manner as the solutions with R m = 0, and so the R[ n = Green's function suffices; 
at late times, once the gas reaches the vicinity of the boundary it exhibits the behavior 
£ oc R~ n (l — ^JR~~JR) in that neighborhood. 



1 The factor aris es from assuming that O is nearly Keplcrian at the radius where the torque g oc d^l/dR = 
(JFrank et al.ll2002h . 
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3.2. Zero Mass Flux at R m > 



We now consider the boundary condition of zero mass flow at a particular radius, 



M(Rin) OC ^ (z/Si? 1 / 2 ) 



R=R- ^ ^ ^ MW'W + 5 (*)^W1} 



0. 



(35) 



y=yu 



From the relations d[x e Ji(x)]/dx = x l Ji_\(x) and d[a^Y^(a;)]/0x = x F£_i(x), we obtain the 
relationship between A and i? corresponding to the boundary condition: 



A(k) 
B(k) 



Je-iM 



(36) 



The solution is then 

POO 

E(R,t) = / C{K)R- n ~ l/i [Ji(kx) Vi_i(«) - Y^kx) J £ _i(«;)]exp 






(37) 

Pringld ( 119911 1 solved the special case n — 1 analytically, and noted the mathematical dif- 
ficulty in deriving a solution for a more general case. We find that the mode weight C(k) 
can in fact be solved for with the use of the recently proved generalized Weber transform 
jzhang fc To^g1l2007h . 



Wi(k,x\ a, 6) 



&e(n) k dn, 



Q 2 e (n;a,b) 

We(n, x; a, b) (f>e(x) x dx. 



M x ) 

The functions Wi(k,x; a, b) and Qj{n; a, b) are defined as follows: 

W((k, x; a, b) = Ji(nx) [aYi(n) + bn Y^(k)] — Y^kx) [aJi(n) + bn J^(ac)] 

= Je(Kx) [(a - ib) Yi(k) + 6k y*_i(«)] - Y £ (kx) [(a - £b) J e (n) + 6k J*_i(«)] 

Q 2 (k; a, 6) = [aY«(«) + 6k ^/(«)] 2 + [oJ^k) + 6k ^(«)] 2 

= [(a - £6) F^(k) + 6k Y^k)} 2 + [(a - lb) J t ( K ) + 6k J^_i(k)] 2 . 



(38) 
(39) 



(40) 



(41) 

Above, J' e and Y' t are the ordinary derivatives of the Bessel functions. If a = 1 and 6 = 0, 
the pair is identical to the ordinary Weber transform (equations 1271 and 1281) . 

The choice a = £ and 6=1 corresponds to the desired boundary condition M(R in , t) = 0. 
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The Green's function is then: 



Gr(il, R , t) 



II 



x 

'0 

x exp 



jR -n-l/4 jR /5/4 jR n-2 
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-2(1-^ 2 - 2 



K 



JIM 

K dn. 



Y?-M) 



(42) 



A specific instance of the above Green's function was derived by iPringld ( 1l99ll ) for the 
case n — 1. We can use equation ( 142]) to reproduce that previous solution by noting that 
J_i/ 2 (x) = — Y"i/ 2 (x) = \Jit /2 x _1 / 2 cosx and y_i/2(x) = Jx/iix) — a/7t/2 x~ 1 / 2 sinx. We 
obtain: 



G(R,R',t) 



TTRin \RJ 

R- 3 / 2 R'R: 

2x72^ 
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-1/2 



^!Ain 
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exp 



cos [k(x — 1)] cos [av(x' — 1)] exp 



k 2 £ 

^j wy.ii 



<1h 



[x — x') t 



+ exp 



(x + x'-2) 2 t 



The only difference between this Green's function and the one for n 
R in (equation [32]) is the sign in between the exponential functions. 



(43) 
1 and zero torque at 



For general values of n, the analytic late-time behavior of equation (142)) turns out to be 
identical to that for the case R[ n = (equation 123]) . This can be confirmed by observing that 
for £ < 1 and small arguments k< 1 and kx <C 1, We(ii,x;£,l) ~ csc(£7r) J_i(kx) J^_i(«) 
and QJ(k;£, 1) ~ csc 2 (£-7r) J 2 _ 1 (k), and therefore the large fraction in equation (T4"2"]) is ap- 
proximately equal to J^(kx) J^(kx'). However, because the disc does not extend to the 
origin for a finite boundary, the integral for the central disc luminosity converges. For the 
^-function initial surface density profile, we obtain: 

t 



L(R<R u t>t v {R)) 



GMM SS 




Rt 



n\- 



t 



2/ t ■ 



(44) 



where again M ss is the mass supply expression defined in §3.1] and R t is the radius where 
t = t u (R), inside which the disc has had sufficient time to approach the asympt otic solution. 
The a bove expression for the disc luminosity is in agreement with the estimate of llvanov et al. 
(119991 ). who considered a zero-flux boundary condition in the context of a thin disc around 
a supermassive black hole binary. 

Figure H] shows the solution for the 5-function initial condition and the zero- flux bound- 
ary condition at R in = Rq/5. The panels showing the mass flow clearly exhibit the desired 
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bo undary conditi on. Note that the case n—\ (panels b and d) is the case solved analytically 
by iPringld (Il99lf ). The n — 1 case, however, leads to a more rapid evolution and steeper 
late-time profiles than solutions with lower values for n; e.g., for the innermost regions of 
circumbinary discs ar ound supermassive black holes , the viscosity is believed to be roughly 
constant with radius (jMilosavljevic &: Phinneyll2005l ; iTanaka &: Menoull2010l ). 



4. Conclusion 

We have presented Green's functions to the equation for viscous diffusion in a thin 
Keplerian accretion disc, in the special case of a power-law viscosity profile v oc R n , for two 
different types of boundary conditions, zero viscous torque or zero mass flow, imposed at 
a finite i nner r adius R m > 0. They are extensions of the elegant analytic solutions derived 
by iLiistl (119521 ) and ILP74I for the same boundary conditions applied at R in = 0. While the 
problem of the finite-radius boundary had been mentioned previously in the literature, to 
the author's knowledge these solutions have not been explicitly pursued, and are presented 
here for the first time. The new solutions can be used to model the time-dependent behavior 
of the innermost regions of accretion discs, where the finite physical size of the central 
objects can significantly affect the observable characteristics of the disc. Whereas the power 
viscously dissipated in the -R; n = solutions diverge, and require manipulation of the profile 
at the disc center to calculate physically plausible disc luminosities, the power for the new 
solutions converge to expressions that are consistent with disc luminosities inferred by other 
(non-Green's function) methods. The solutions presented here complement the numerous 
approximate solutions and numerical treatments in the literature^ 

The integral transforms used to derive the solutions are applicable to a wide class of 
boundary conditions, and may be applicable to astrophysical thin-disc systems and con- 
figurations not considered here. Because the generalized Weber transform by its nature is 
applicable to many second-order differential equations with intrinsic cylindrical symmetry, 
they may also prove to be useful in solving other mathematical equations in astrophysics 
and other fields. 



2 For example, ICannizzo et al.l (|l990l ) studied the accretion of a tidally disrupted star onto a black hole 
via numerical solutions and analytic self-similar solutions. The problem of a thin disc with M = at a finite 



radius was discussed for the non- linear case v oc £"V™ by lPringld ( 19911 ) and llvanov et al.l < 19991 ). with both 
papers providing numerical solutions as well as analytic approximations. 
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n = 0.1 



g(0,t) = 



n=l 



w 

J-J 



•2 




Fig. 1.— The solution E(i£, £), from ILP74I . and the corresponding radial mass inflow rate 
M(R,t), for a zero-torque boundary condition imposed at R — and a ^-function initial 
profile T>(R,t = 0) = Eq-Ro'K-R — -Ro) ; where the quantities E and R are arbitrary. The 
viscosity is a radial power law with v oc i? n . Panels on the left side (a and c) show solutions 
for n = 0.1, and those on the right (b and d) show solutions for n = 1. Values for t are 
in units of the viscous time at Rq, t v $ = (2/3)Rq/u(R ). We have normalized M to the 
quantity M = 37ri/(i?o)S . At late times, the solution has the behavior E oc R~ n and the 
mass-flow profile M becomes flat near the origin. 
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n = 0.1 



M(0,t) = 



n=l 



w 
\ 

J-J 



•5§ 




Fig. 2.— Same as Figure (TJ except that the boundary condition is M = at R = 0. 
Again, the scales S and Rq are arbitrary. Whereas in the zero-torque case the total mass 
in the disc monotonically decreases due to mass loss at the origin (onto the black hole or 
star), the solutions in this figure conserve mass. At late times, the solution has the behavior 
S oc i?"™" 1 / 2 . Gas initially piles up near the origin because of the boundary condition before 
gradually spreading outward; note that E at inner radii decreases from t = 0.3£„ i0 to t = t u> Q. 
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Fig. 3. — Same as as Figure [U except that the zero-torque boundary condition is applied at 
a finite radius -Rj n = Rq/5. As gas flows n ear the inner boundary, it exhibits the well-known 



behavior E oc R~ n (l - ^R in /R) of [LP74. 
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Fig. 4. — Same as as Figure [21 except that the zero-flux boundary conditio n is appl i ed at a 
finite radius R[ n = Rq/5. Note that the n = 1 case was solved analytically by iPringld ( 1l99ll ). 



